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The  main  thrust  of  our  research  program  has  been  the  development  and 
applications  of  asymptotic  and  perturbation  for  analyzing  the  stability  and 
dynamics  of  elastic  structures  and  fluid  flows.  The  work  is  summarized  in 
the  following  papers  which  have  been  published,  accepted  for  publication, 
submitted  for  publication,  or  are  in  preparation  for  publication. 

1.  Sinay,  L.R.,  and  Reiss,  E.L.,  Perturbed  Panel  Flutter:  A  Simple  Model, 

AIAA  J.,  19  (1981),  pp.  1476-1483. 

The  effects  of  small  periodic  disturbances  on  the  reponse  of  a  two- 
degree-of- freedom,  non-conservative,  mechanical  system  are  analyzed.  The 
system  is  a  simple  model  for  panel  flutter.  The  distrubance  simulates  the 
pressure  fluctuations  of  a  turbulent  boundary  on  the  panel.  Asymptotic 
expansions  of  the  solutions  are  obtained  for  small  amplitudes  of  the  dis¬ 
turbance.  The  qualitative  features  of  the  response  depend  on  the  prescribed 
variation  of  the  frequency  of  the  disturbance  with  the  magnitude  of  the 
non-conservative  applied  force.  The  disturbance  can  induce  a  smooth 
transition  to  the  fluttering  states  of  the  rods;  or  it  may  induce  jump 
transitions.  The  results  suggest  a  possible  technique  for  delaying  panel 
flutter,  by  imposing  a  periodic  forcing  function  with  an  appropriate 
frequency . 

2.  Watson,  J.  G.,  and  Reiss,  E.  L.,  A  Statistical  Theory  of  Imperfect 
Bifurcation,  SIAM  J.  Appl.  Math.,~42  (1982),  pp.  135-148. 

An  "honest"  statistical  method  is  presented  to  analyze  the  effects  of 
imperfections,  and  other  disturbances  on  the  bifurcation  of  solutions  of 
nonlinear  problems.  First,  uniformly  valid  asymptotic  approximations  of 
the  solutions  are  obtained  for  any  realization  of  the  imperfections. 

The  approximations  are  valid  as  the  magnitude  of  the  imperfections 
approaches  zero.  The  statistical  properties  of  the  solutions  are  then 
deduced  directly  from  these  approximations,  for  specified  statistics  of 
the  imperfections.  For  simplicity  of  presentation,  the  imperfections 
are  taken  as  zero-mean,  wide-sense  stationary,  Gaussian  random  processes. 

The  statistical  analysis  is  elementary.  It  provides  easily  analyzed 
results  for  the  expected  values  and  variances  of  the  solutions.  Confidence 
limits  are  also  given. 

3.  Dellwo,  D.,  Keller,  H.  B.,  Matkowsky,  B.J.,  and  Reiss,  E.  L.,  On  the 
Birth  of  Isolas ,  SIAM  J.  Appl.  Math.,  42  (1982),  pp.  956-963. 

Isolas  are  isolated,  closed  curves  of  solution  branches  of  nonlinear 
problems.  They  have  been  observed  to  occur  in  the  buckling  of  elastic 
shells,  the  equilibrium  states  of  chemical  reactors,  and  other  problems. 

In  this  paper,  we  present  a  theory  to  analytically  describe  the  structure 
of  a  class  of  isolas.  Specifically,  we  consider  isolas  that  shrink  to  a 
point  as  a  parameter  k  of  the  problem,  approaches  a  critical  value  Kq. 

The  point  is  referred  to  as  an  isola  center.  Equations  that  characterize 
the  isola  centers  are  given.  Then  solutions  are  constructed  in  a  neighbor¬ 
hood  of  the  isola  centers,  by  a  perturbation  expansion  in  a  small  parameter 
c,  that  is  proportional  to  (k  -  Kgfy,  with  ^appropriately  chosen.  The 
theory  is  applied  to  problems  in  chemical  reactor  theory,  and  nonlinear 
oscillations. 
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4*  Grotberg,  J.B.,  and  Reiss,  E.  L.,  A  Subsonic  Flutter  Anomaly,  J.  Sound 
and  Vibration,  80  (1982),  pp.  444-446. 

Experiments  on  the  flutter  of  cylindrical  tubes  conveying  subsonic  flows 
reveal  that  the  tubes  first  partially  collapse  into  a  static;  nearly 
flat  state  at  a  critical  flow  velocity.  Then  at  a  second  and  larger 
critical  flow  velocity  the  tube  flutters  about  this  flat  state.  These 
results  are  at  variance  with  existing  theoretical  studies  which  show  that 
the  nearly  flat  state  loses  stability  by  divergence.  In  these  theories 
the  fluid  is  assuaed  to  be  invlscid  and  incompressible.  He  have  resolved 
this  anoaaly  by  including  fluid  damping  in  the  model.  Then  the  nearly 
flat  state  loses  stability  by  flutter.  It  is  low  frequency  flutter  for 
small  damping. 


5.  Strumolo,  G«,  and  Reiss,  E.  L.,  Poiseuille  Channel  Flow  with  Driven  Walls, 
J.  Fluid  Mech.,  submitted. 

The  effects  of  a  prescribed  wall  antion  on  the  nonlinear  stability  of 
Poiseuille  channel  flow  are  studied  by  an  asymptotic  method.  The  motion 
represents  a  traveling  wave  in  the  upper  wall  of  the  channel.  It  can 
be  considered  either  as  a  disturbance  to  the  flow  that  results  from 
experimental  imperfections,  or  as  an  externally  imposed  motion.  The 
frequency  of  this  disturbance  depends  on  the  Reynolds  number  of  the  flow. 
In  the  classical  Poiseuille  channel  flow  problem,  the  walls  are  assumed 
to  be  rigid.  Then  a  periodic  solution  bifurcates  from  the  laminar, 
Poiseuille  flow  at  the  critical  Reynolds  number,  R  .  In  the  resonance 
case,  the  wall  motion  destroys  the  bifurcation.  Tne  transition  from 
the  laminar  state  then  occurs  by  jumping  at  new  critical  Reynolds  numbers. 
These  Reynolds  numbers  either  exceed,  or  are  below  R  ,  depending  on  the 
variation  of  the  wall  motion  frequency  with  the  Reynolds  number.  Thus 
the  wall  motions  can  stabilise  or  destabilise  the  laminar  flow,  and  hence 
they  can  be  used  to  control  the  transition  to  turbulence.  In  the  non¬ 
resonance  cases,  the  bifurcation  is  preserved  and  the  critical  Reynolds 
number  is  slightly  perturbed. 


6.  Reiss,  E.  L. ,  Cascading  Bifurcation.  SIAM  J.  Appl.  Math.,  43  ( 1983), pp. 57-65. 

Sequential  bifurcation  of  solutions  of  nonlinear  equations,  as  the  bifur¬ 
cation  parameter  increases,  is  called  cascading  bifurcation.  It  has  been 
proposed  as  a  mechanism  to  describe  the  transition  from  laminar  to 
turbulent  fluid  flows,  and  as  a  mechanism  for  chemical  and  biochemical 
morphogenesis  and  pattern  formation.  The  creation  of  cascading  bifurcation 
by  the  splitting  of  multiple  primary  bifurcation  points  is  described. 

A  perturbation  method  is  employed  in  the  analysis.  Primary,  secondary, 
tertiary,  etc.  bifurcation  points  and  states  are  determined  for  a 
specific  nonlinear  boundary  value  problem. 
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7.  Slnay,  L.  R. ,  and  Reiss,  E.  L. ,  Secondary  Transltlona  in  Panel  Flutter: 

A  Single  Model,  to  be  submitted. 

A  simple,  two  degree  of  freedom  mechanical  model  of  panel  flutter  is 
presented.  A  perturbation  method  Is  employed  to  determine  the  secondary 
bifurcation  of  flutter  states  from  divergence  states,  and  divergence 
states  from  flutter  states.  The  latter  suggests  a  new  method  for  controlling 
flutter  by  allowing  small  amplitude  flutter  and  then  increasing  the  flow 
velocity  until  secondary  bifurcation  into  a  divergence  state  occurs. 

8.  Grotberg,  J.  B.,  and  Reiss,  E.  L. ,  Subsonic  Flapping  Flutter.  J.  Sound 
and  Vibration,  submitted. 

A  mathematical  model  of  two-dimensional  flow  through  a  flexible  channel 
is  analyzed  for  its  stability  characteristics.  Linear  theory  shows 
that  fluid  viscosity,  modelled  by  a  Darcy  friction  factor.  Induces 
flutter  instability  when  the  dimensionless  fluid  speed,  S,  attains  a 
critical  flutter  speed,  SD.  This  is  in  qualitative  agreement  with 
experimental  results,  and  it  is  at  variance  with  previous  analytical 
studies  where  fluid  viscosity  was  neglected  and  divergence  instability 
waa  predicted.  The  critical  flutter  speed  and  the  associated  critical 
flutter  frequency  depend  on  three  other  dimensionless  parameters:  the 
ratio  of  fluid  to  wall  damping;  the  ratio  of  wall  to  fluid  mass;  and 
the  ratio  of  wall  bending  resistance  to  elastance.  Non-linear  theory 
predicts  stable,  finite  amplitude  flutter  for  S  >  SQ,  which  Increases 
in  frequency  and  amplitude  as  S  increases.  Both  symmetric  and  antisymmetric 
modes  of  deformation  are  discussed. 


9.  Magnan,  J.,  and  Reiss,  E.  L.,  Double-Diffusive  Convection  and  1 -Bifurcation, 
to  be  submitted. 

We  consider  convection  in  a  rectangular  box  where  two  "substances"  such 
as  temperature  and  a  solute  are  diffusing.  The  solutions  of  the 
Bousslnesq  theory  depend  on  the  thermal  and  solute  Rayleigh  numbers  Rj 
and  R8  in  addition  to  other  geometrical  and  other  fluid  parameters.  The 
conduction  state  is  unstable  with  respect  to  steady  (periodic)  convection 
states  if  Rg  is  sufficiently  small  (large).  The  boundary  between  steady 
and  periodic  convection  occurs  at  a  critical  value  Rs  ■  Rs.  The  linearized 
theory  at  Rg  »  Rg  is  characterized  by  the  frequency  «  ■  0  appearing  as  a 
root  of  algebraic  multiplicity  two  and  geometrical  multiplicity  one._ 
Asymptotic  approximations  of  the  solutions  are  obtained  for  Rg  near  R_ 
by  the  Poincare-Linatedt  method.  It  la  found  that  a  periodic  (steady) 
solution  bifurcates  supercritlcally  (subcrltiCally)  from  the  conduction 
state  at  Ry  •  R£(Rf),  where  Ry  <  R|.  The  periodic  branch  joins  the  steady 
branch  with  an  "infinite  period  bifurcation"  at  Ry  *  Rb,  where  Ry  <  R.  <  Ry'. 
The  shape  of  the  resulting  bifurcation  diagram  suggests  the  term,  bifurcation. 
The  infinite  periodic  bifurcation  corresponds  to  a  heteroclinic  orbit  in 
the  appropriate  amplitude  phase  plane.  The  periodic  (steady)  convection 
states  are  stable  (unstable),  aa  we  demonstrate  by  solving  the  initial 
value  problem  employing  the  nultl-scale  method. 


-4 


10.  Erneux,  T. ,  and  Reiss,  E.  L. ,  Brussellator  Isolas.  SIAN  J.  Appl.  Math., 
in  press. 

The  Brussellator  is  a  simple  chemical  model  describing  pattern  formation 
by  bifurcation  of  solutions.  For  a  one-dimensional  system,  the  bifurcation 
parameter  is  related  to  the  ratio  of  the  square  of  the  size  of  the  system 
to  a  diffusion  coefficient.  It  has  been  observed  from  numerical  computations, 
that  there  are  closed  branches  of  steady  state  solutions,  which  are  called 
Isolas,  that  connect  neighboring  bifurcation  points.  In  addition,  these 
isolas  depend  on  a  parameter  B.  As  B  approaches  a  critical  value  B°  the 
neighboring  bifurcation  points  coalesce,  so  that  the  isola  shrinks  to  a  point. 
We  employ  a  perturbation  method  to  obtain  asymptotic  expansions  of  the  Isolas 
for  3  near  B°.  Implications  of  the  results  for  pattern  formation  are 
discussed. 


11.  Erneux,  T.,  and  Reiss,  E.  L. ,  Splitting  of  Steady  Multiple  Eigenvalues 
May  Lead  to  Periodic  Bifurcation.  SIAM  J.  Appl.  Math.,  accepted. 

A  general  bifurcation  problem  is  considered  that  depends  on  two  parameters 
In  addition  to  the  bifurcation  parameter  X.  It  is  assumed  that  all  primary 
bifurcation  states  correspond  to  steady  solutions  and  that  they  branch 
supercritlcally.  Then  it  is  shown  that  for  a  range  of  system  parameters, 
and  near  a  triple  primary  bifurcation  point  the  following  cascade  of  bifur¬ 
cations  from  the  minimum  primary  bifurcation  state  is  possible.  As  \ 
increases  there  is  secondary  and  then  tertiary  bifurcation  to  steady  states, 
and  finally  Hopf  bifurcation  at  a  quarternary  bifurcation  point.  Related 
transitions  have  been  observed  experimentally  in  thermal  convection  and 
other  hydrodynamic  stability  problems.  In  addition,  we  show  that  Hopf 
bifurcation  near  a  double  primary  bifurcation  point  is  not  possible  when 
both  primary  states  near  the  double  point  bifurcate  supercritlcally. 

However,  it  is  possible  near  such  a  double  bifurcation  point  if  Imperfections 
are  Included  in  the  formulation,  as  we  demonstrate. 


12.  Kriegsmann,  G.  A.,  Radiation  Conditions  for  Wave  Guide  Problems.  SIAM 
Scl.  Stat.  Comput.,  3  (1982),  pp.  318-326. 

An  incident  mode  propagates  down  a  two-dimensional  wave  guide  until  it 
strikes  a  localized  obstruction  which  creates  reflected  and  transmitted 
waves.  The  numerical  determination  of  these  waves  is  difficult  because 
the  classical  radiation  condition  does  not  apply  for  an  infinite  wave 
guide.  In  this  note  we  derive  a  sequence  of  "localized  radiation 
conditions"  which  can  be  applied  a  few  wavelengths  away  from  the  scat¬ 
tering  object.  These  conditions  allow  us  to  numerically  solve  the 
Helmholtz  equation  on  a  finite  domain. 
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13.  Krlegs— nn ,  C.  A.,  Exploiting  the  Limiting  Aaglitudg  Principle  r> 

Hu— rlcallv  Solve  Scattering  Problems.  Wave  Motion  4  (1982),  pp.  371-380. 

A  nu— rical  method  for  solving  reduced  wave  equation*  1*  presented. 

The  technique  1*  basically  a  relaxation  ache—  which  exploit*  the 
Halting  amplitude  principle.  A  modified  radiation  condition  at  "infinity" 
1*  also  given.  The  method  is  tested  on  two  model  proble— :  the  scattering 
of  plane  shallow  water  waves  off  shoals  and  the  scattering  of  plane 
acoustic  waves  off  a  sound-soft  cylinder  imbedded  between  two  homogeneous 
but  different  half  spaces.  The  nu— rical  solutions  exhibit  correct 
refractive  and  diffractive  effects  at  —derate  frequencies. 


14.  Magnan,  J.  F.,  and  Goldstein,  R.  A.,  Perturbed  Bifurcation  of  Stationary 
Strlatlons  in  a  Contaminated.  Nonunlform  Plasma.  SIAM  J.  Appl.  Math., 

43  (1983),  pp.  16-30. 

A  cylindrical,  weakly  ionized  and  collision  '  minated  neon  plasma  can 
be  described  by  a  system  of  nonlinear,  paraboxic  reaction-diffusion 
equations  for  the  electron  and  metastable  atom  axial  densities.  The 
equations  exhibit  a  bifurcation  from  a  uniform  to  a  striated  state  at 
a  critical  length  of  the  plasma  column.  The  sharp  transition  between 
states  predicted  by  the  theory  is  in  contrast  with  the  smooth  transition 
observed  in  experiments.  He  apply  the  theory  of  singular  perturbation 
of  bifurcations  to  show  that  s— 11  inhomogeneities  in  the  plasma ,  such 
as  those  caused  by  nonuniform  heating  and  contamination,  are  sufficient 
to  qualitatively  explain  the  experl— ntal  results.  We  observe  that  a 
steady,  axial  — gnetlc  field  in  the  plasma  can  also  produce  a  smooth 
transition. 


15.  Magnan,  J.  F.,  and  Huerta,  M.  A.,  Spatial  Structures  in  Plasma s  with 

Metastable  States  as  Bifurcation  Pheno— na.  Physical  Review  A,  26  (1982), 
pp.  539-555. 

The  for— tion  of  orderly,  steady,  spatial  structures  in  rf-heated  plasmas 
that  have  metastable  states  is  studied  as  the  development  of  an  instability 
of  the  solution  to  the  reaction-diffusion  equations  for  the  plas—  species. 
Plas— -species  (electrons,  ions,  and  — tastable  ato— )  source  functions 
are  given  in  ter—  of  the  electron  energy  distribution  in  strong  electric 
fields  and  of  the  various  collision  cross  sections.  The  — thods  of  two- 
tl—  expansions  and  bifurcation  theory  are  used  to  calculate  the  non- 
uniform  states  that  bifurcate  from  the  unstable  uniform  steady  state 
for  straight  and  toroidal  tubes.  The  sinusoidal  density  variations  give 
the  discharge  the  appearance  of  glowing  balls  of  gas.  This  pattern 
for— tion  is  analogous  to  the  morphogenesis  found  in  certain  chemical 
reactions  and  biological  processes. 


16.  Grotberg,  J.  B.,  Louhi,  B. ,  and  Reiss,  E.  L. ,  Secondary  Bifurcation  of 
Quaslperlodic  Solutions  Can  Lead  to  Period  Multiplication,  submitted. 

It  la  shown  using  perturbation  and  asymptotic  methods  that  secondary 
bifurcation  of  quasl-periodlc  solutions  from  periodic  solutions  occurs 
for  a  model  problem.  The  model  Is  a  coupled  system  of  two  van  der  Pol- 
Duffing  oscillators.  For  special  values  of  the  detuning  parameters 
the  secondary  states  are  periodic.  Then  periodic  multiplication  of 
solutions  can  occur  at  the  secondary  bifurcation  point. 


17.  Krlegsmann,  G.  A.,  and  Reiss,  E.  L.,  Low  Frequency  Scattering  by  Local 
Inhomogeneities ,  SIAM  J.  Appl.  Math.,  4_3  (1983). 

An  asymptotic  expansion  which  is  uniformly  valid  in  space  is  obtained 
for  the  low  frequency  scattering  of  a  plane  wave  incident  on  a  localized 
inhomogen«lty.  The  scattering  region,  which  may  be  simply  or  multiply 
(collection  of  scatterers)connected,  has  a  characteristic  length  which 
is  small  compared  with  the  wave  length  of  the  incident  wave.  The 
index  of  refraction  n  is  unity  outside  the  scattering  region  and  it  is 
arbitrary  Inside  the  region.  The  method  of  matched  asymptotic  expansions 
is  used  in  the  analysis.  The  Born  approximation  is  shown  to  agree  with 
the  uniform  expansion  in  the  far  and  the  near  fields.  The  leading  term 
in  the  uniform  expansion  is  a  linear  functional  of  1  -  n2.  Thus,  statistics 
of  the  scattered  field  are  easily  evaluated  from  the  statistics  of  n, 
when  n  is  a  random  process.  The  method  is  then  applied  to  the  low  fre¬ 
quency  scattering  of  an  acoustic  plane  wave  by  localized  inhomogeneltles 
in  the  density  and  the  index  of  refraction.  Finally,  the  scattering  by 
a  plane  which  is  acoustically  hard  except  for  a  small  impedence  spot, 
is  analyzed  by  the  same  method. 
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The  stability  of  the  null  state  for  a  nonlinear  Burgers  system  is  examined. 
The  results  Include  (1)  an  energy  estimate  for  global  stability  for  states 
involving  arbitrary  modulation  in  time,  and  (ii)  an  analysis  of  the  bifur¬ 
cation  from  the  null  state  for  slow  modulations.  For  the  slow  modulations 
it  is  determined  that  the  amplitude  A(t)  of  the  bifurcated  disturbance 
velocity  satisfies  a  Landau-type  equation  with  time-dependent  growth 
rate  0(t).  Particular  attention  is  given  to  periodic  and  quaslperlodic 
modulations  of  the  system,  which  lead  to  analogous  behavior  in  6(t).  For 
each  of  these  oscillatory-type  modulations,  it  is  found  that  A^(t)  has 
the  same  long-time  mean  value  as  the  unmodulated  case,  implying  no 
alteration  of  the  final  mean  kinetic  energy.  Applications  to  various 
fluid-dynamical  phenomena  are  discussed. 
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A  method  Is  developed  for  solving  an  inverse  problem  for  the  wave  equation 
with  potential  where  the  object  is  to  find  the  potential  given  Cauchy 
data  on  a  time- like  surface.  The  computation  is  carried  out  with  one 
space  variable  by  an  iterative  procedure.  The  point  of  this  method  is 
that  it  can  be  extended  to  higher  dimensions  in  principle/  A  coarse 
mesh  finite  difference  scheme  is  used  which  yields  fair  accuracy. 
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Flows  and  the  Formation  of  Caustics,  submitted. 

The  propagation  of  acoustic  waves  from  a  high  frequency  line  source  in 
a  two  dimensional  parallel  shear  flow  adjacent  to  a  rigid  wall  is 
analysed  by  a  ray  method.  The  leading  term  in  the  resulting  expansion 
is  equivalent  to  the  geometrical  acoustics  theory  of  classical  wave 
propagation.  It  is  shown  that  energy  from  the  source  is  radiated  either 
directly  to  the  far  field,  or  by  first  reflecting  from  the  wall.  In 
addition,  energy  is  trapped  in  a  channel  adjacent  to  the  wall  and  down¬ 
stream  from  the  source.  The  rays  in  this  channel  form  an  infinite 
sequence  of  caustics  progressing  downstream.  Since  the  geometrical 
acoustics  approximation  is  invalid  on  and  near  caustics,  a  boundary 
layer  method  is  employed  to  determine  the  acoustics  field  near  the 
caustics.  It  is  shown  that  the  amplitude  of  the  field  on  and  near  the 
caustics  is  k2'3  larger  than  the  geometrical  acoustics  field  for  large  k. 

Here  K  is  a  dimensionless  wave  number  of  the  source.  Moreover,  the  vorticity 
of  the  acoustics  field  in  the  caustic  regions  is  k''^  larger  than  the 
geometrical  acoustics  field.  The  possible  significance  of  these  results 
for  vehicle  self-noise  and  the  formation  of  turbulent  spots  in  the  sub-layer 
of  a  turbulent  boundary  layer  is  discussed. 


